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Abstract 

The inert doublet model, a minimal extension of the Standard Model by a second Higgs 
doublet, is one of the simplest and most attractive scenarios that can explain the dark matter. 
In this paper, we demonstrate the existence of a new viable region of the inert doublet 
model featuring dark matter masses between Mw and about 160 GeV. Along this previously 
overlooked region of the parameter space, the correct relic density is obtained thanks to 
cancellations between different diagrams contributing to dark matter annihilation into gauge 
bosons (W + W~ and Z°Z°). First, we explain how these cancellations come about and show 
several examples illustrating the effect of the parameters of the model on the cancellations 
themselves and on the predicted relic density. Then, we perform a full scan of the new viable 
region and analyze it in detail by projecting it onto several two-dimensional planes. Finally, 
the prospects for the direct and the indirect detection of inert Higgs dark matter within this 
new viable region are studied. We find that present direct detection bounds already rule out 
a fraction of the new parameter space and that future direct detection experiments, such as 
XenonlOO, will easily probe the remaining part in its entirety. 

1 Introduction 

The identification of the dark matter particle is one of the most challenging problems in astropar- 
ticle physics today. We know that dark matter exists, for the evidence in its favor is varied and 
overwhelming. It comes, among others, from rotation curves in galaxies, gravitational lensing 
effects in clusters of galaxies, the large scale structure of the Universe, and the Cosmic Microwave 
Background (CMB). The anisotropics in the CMB, in fact, indicate that while dark matter ac- 
counts for about 23% of the energy density of the Universe, baryonic matter reaches only around 
4% pp. Most of the matter in the Universe, therefore, is in the form of dark matter. And yet we 
still ignore what dark matter consists of. 

Currently there is a huge experimental effort aimed at detecting the dark matter particle and 
determining its main properties. Dark matter particles can be detected, for example, via nuclear 
recoil as they scatter off nuclei in direct detection experiments. In recent years, direct detection 
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experiments such as CDMS [2] or Xenon [3] have been continuously improving the constraints on 
the spin-independent dark matter-nucleon cross section, and they have started to probe interesting 
regions of the parameter space of common dark matter models. The indirect detection of dark 
matter, which relies on the annihilation of dark matter particles rather than on their scatterings, is 
another promising avenue toward the identification of the dark matter. Dark matter annihilations 
may produce photons, antimatter (positrons and antiprotons mainly) and neutrinos that can be 
detected as exotic components in cosmic rays. At present, several experiments including FERMI 
(photons) [1], PAMELA (antimatter) [5] and ICECUBE (neutrinos) [5] are already taking data 
and setting some constraints on the dark matter properties. Finally, dark matter particles or other 
particles associated with it may also be produced and detected in colliders such as the LHC. It 
may be that the LHC soon finds evidence of physics beyond the Standard Model that sheds some 
light on the dark matter puzzle. 

In parallel with this experimental effort, there has been a lot of interest in the theoretical 
aspects of dark matter. Even though extremely successful in describing collider experiments, the 
Standard Model of particle physics does not include any viable dark matter candidate. Dark 
matter, therefore, provides compelling evidence for the existence of new physics, of physics be- 
yond the Standard Model. Over the years, many dark matter candidates have been proposed in 
different scenarios for physics beyond the Standard Model. They include the neutralino [7] and 
the gravitino [8] in supersymmetric models, the lightest Kaluza-Klein particle in Universal Extra 
Dimensional models [5], as well as the singlet scalar [TU] or the inert Higgs [TTHT5] that appear 
in minimal extensions of the Standard Model. Thanks to the simplicity of the underlying model 
and to the rich phenomenology associated with it, the inert Higgs -the lightest odd particle of the 
inert doublet model- has earned a special place among them as one of the most attractive dark 
matter candidates considered in the recent literature. 

In the inert doublet model, a Higgs doublet H2, odd under a new Z2 symmetry, is added to 
the Standard Model particle content. The lightest inert (odd) particle, H°, turns out to be stable 
and hence a suitable dark matter candidate. After being introduced in [TT], this model has been 
extensively studied in a number of recent works [P2TT18] . In [12] , it was extended with odd right- 
handed neutrinos so as to explain neutrino masses. Then, in [I5j . it was proposed as a possible 
way to increase the Higgs mass without contradicting electroweak precision data. Regarding dark 
matter, the main features of the model were systematically analyzed in [14] . Later on, in |15j . 
it was shown that the annihilation of inert Higgs dark matter can give rise to significant gamma 
ray lines. Recently, it was demonstrated that, for masses between 60 GeV and M\y, dark matter 
annihilations are usually dominated by three-body final states consisting of a real and a virtual 
W [ElEO], rather than by the two-body final states considered in previous works. A result that 
has important implications for the determination of the viable parameter space of the model and 
for the detection prospects of inert Higgs dark matter. 

These various studies have established that in the inert doublet model the dark matter con- 
straint can only be satisfied within specific regions of the parameter space. Three viable regions 
have been identified: a low mass regime, with M H o < 50 GeV, where dark matter annihilates 
exclusively into fermion pairs; an intermediate mass regime, with M H o < My/, where dark matter 
annihilations are dominated by the final states WW* and bb, and where coannihilations with 
A may also play a role in the determination of the relic density; and a heavy mass regime, 
Mho > 500 GeV, where the annihilation into gauge bosons is dominant and the mass splitting 
between the inert particles is small. 

We demonstrate in this paper the existence of another viable region of the inert doublet model, 
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a region that had been overlooked in previous analysis. This region features dark matter masses 
between My/ and about 160 GeV and is compatible with the relic density constraint thanks to 
cancellations among the different diagrams that contribute to dark matter annihilation into gauge 
bosons. It is our goal to introduce this new viable region and to analyze its main properties. We 
determine its position in the multidimensional parameter space of the inert doublet model and 
study its implications for the direct and indirect detection of inert Higgs dark matter. Remarkably, 
we find that this new viable region will be entirely and easily probed by the next generation of 
direct dark matter detection experiments such as XenonlOO. 

In the next section, the inert doublet model is briefly reviewed, with particular emphasis on 
the different bounds that must be taken into account. Then, we describe in some detail the known 
viable regions of the model, those that are compatible with the observed dark matter abundance 
and with all other experimental and theoretical constraints. In section|4l we illustrate, analytically 
and numerically, the existence of cancellations in the annihilation of inert Higgs dark matter into 
W + W~ and Z°Z° that allow to satisfy the relic density constraint for dark matter masses above 
My/. The effect of the different parameters of the model on these cancellations and on the 
predicted relic density are studied, and several examples of models within the new viable region 
are provided. In section [SJ we use Markov Chain Monte Carlo techniques to scan the parameter 
space of the inert doublet model and obtain a large sample of models within the new viable region. 
We then determine the main features of this region by projecting this multidimensional cloud of 
viable models onto several two-dimensional planes. The issue of the fine-tuning required to obtain 
the correct relic density is also addressed, and it is found that these models are not particularly 
fine-tuned. Finally, we consider, in section [51 the detection prospects of inert Higgs dark matter 
in direct and indirect detection experiments. 



2 The inert doublet model 

The inert doublet model is a simple extension of the Standard Model with one additional Higgs 
doublet i?2 and an unbroken Zi symmetry under which if 2 is odd while all other fields are even. 
This discrete symmetry prevents the direct coupling of H2 to fermions and, crucial for dark matter, 
guarantees the stability of the lightest inert particle. The scalar potential of this model is given 

by 
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+ \i\H\H2\ 2 + ^ [(H\H 2 ) 2 + h.c] , (1) 

where Hi is the Standard Model Higgs doublet, and Xi and /i^ are real parameters. As in the 
Standard Model, Ai and u? can be traded for the Higgs mass, Mh, and the vacuum expectation 
value of Hi, v = 246 GeM±|. Four new physical states are obtained in this model: two charged 
states, H^, and two neutral ones, H° and A . Either of them could be dark matter. In the 
following we choose H° as the lightest inert particle, M^ < M^ ,M^ ± , and, consequently, as 
the dark matter candidate. After electroweak symmetry breaking, the inert scalar masses take 
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the following form 

M 2 H± =n\ + ±\ 3 v 2 , (2) 

M| =M2 + ^3 + A 4 + A 5 )w 2 , (3) 

M 2 A o =M2 + ^(A 3 + A 4 -A 5 )t; 2 . (4) 

Notice that A5 determines the mass splitting between Mjj° and M A o and is necessarily negative as 
M A a > Mho. The mass splitting between the charged and the neutral components, on the other 
hand, is due to both A4 and A5. It is convenient to trade for /i 2 the parameter Ai defined by 

, _ A3 + A4 + A 5 _ M H o — fi 2 

<\l = = — 5 ■ ^ 

2 tr 

As we will see, this parameter, which determines the interaction term between a pair of H and 
the Higgs boson, will play a crucial role in the cancellations discussed in Section |H The free 
parameters of the inert doublet model can then be taken to be the following: 

M H , M H ± , M A o , M h , and X L . (6) 

From them the Aj and [i 2 parameters appearing in the scalar potential, equation ((TJ), can be 
reconstructed via equations (J2J5)). 

These parameters are subject to a number of theoretical and experimental constraints -see [13] 
and 14 . To prevent the breakdown of perturbation theory, the \ parameters cannot be very 
large. We demand accordingly that 

|A,| < 4tt . (7) 

The requirement of vacuum stability, in addition, imposes the following conditions on the \ 
parameters: 

Ai,A a >0, A 3 ,A 3 + A 4 - |A 5 | > -2s/MM- (8) 

Since A2 only affects the quartic interaction term among the inert particles, the phenomenology 
of the model is rather insensitive to its precise value. For definiteness we set A2 = 3.0 throughout 
this work and use the above equation as a constraint on the remaining parameters of the model. 

Accelerator bounds should also be taken into account. The mass of the charged scalar, M H ±, 
is constrained to be larger than 70 — 90 GeV [21] ; M H o + M A o must exceed Mz to be compatible 
with the Z°-width measurements; and some regions in the plane (Mjjo, M40) are constrained by 
LEP II data [T7] . Another important restriction on the parameter space of the inert doublet 
model comes from electroweak precision data. The inert doublet, If 2, contributes to electroweak 
precision parameters such as S and T. It turns out that its contribution to S is always small [13 , 
so it can be neglected. Its contribution to T, AT, was computed in [13] and is given by 

AT = Ta 2 2 i F i M H± , M A o ) + F(M H ± , M H o ) - F(M A o , M H o )] (9) 



where 
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This contribution originates in the A^s terms in the potential, for those are the terms breaking 
the custodial symmetry. The Higgs mass, Mh, also affects T via [22] 

Th ~ ^In^- (ID 

8vrcos 2 6W M z 

For large Mh this contribution violates current experimental constraints. In fact, if no new physics 
affects the electroweak precision parameters, Mh should be smaller than 186 GeV at 95%C.L. In 
the inert doublet model, however, there is certainly physics beyond the Standard Model, so it is 
possible to have a heavier Higgs boson provided that the inert particles give rise to a positive 
compensating AT. In fact, the Higgs boson can be as heavy as 600 GeV in the inert doublet 
model [T5] , 

It is useful to differentiate a light and a heavy Higgs regime in this model. In the heavy Higgs 
regime, Mh > 200 GeV, AT should be positive and equation © implies that M H ± > M A o. In 
the light Higgs regime (Mh < 200 GeV), on the contrary, there is no correlation between M H ± 
and Mao. 

In our numerical analysis we implement the electroweak precision constraint by requiring that 

- 0.1 < AT + T h < 0.2 , (12) 

in agreement with present bounds |23j . 

Finally, since dark matter is the most compelling motivation of the inert doublet model, we 
also constrain the relic density of inert Higgs dark matter to be compatible with the observed dark 
matter density. Explicitly we impose that 

0.09 < Mi 2 < 0.13, (13) 

which is consistent with the WMAP measurement pQ. To ensure that the relic density of in- 
ert Higgs dark matter is accurately computed, we make use of the micrOMEGAs package [53] 
throughout this work. 

As expected, all these constraints can be satisfied only within certain regions of the parameter 
space, the so-called viable regions. Before describing the new viable region of the inert doublet 
model -the main topic of this paper-, we will have a look in the next section to the previously 
known viable regions. 



3 The known viable regions 

In previous works [T3JIII1I2S1I2S] > ^ has been found that the dark matter constraint can only be 
satisfied for restricted values of Mfjo. Three viable regions can be distinguished: a small mass 
regime, with Mjjo < 50 GeV; a large mass regime, where Mjjo > 500 GeV; and an intermediate 
mass regime where 50 GeV > M H o < Mw 

In the small mass range, A and must be decoupled from H° in order to satisfy collider 
and precision test constraints. Thus, the inert doublet model becomes equivalent (for dark matter 
purposes) to the singlet scalar model of dark matter [TU] . Dark matter annihilations proceed in 
this case through Higgs mediated diagrams into light fermions and are therefore dominated by the 
bb final state. Interesting phenomenological issues related to this mass range, including prospects 
for the direct and the indirect detection of dark matter, were studied in [2"BT - |3"0"] . Recently, this 
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small mass regime of the inert doublet model, in particular the region M H o ~ 8 GeV, received 
some attention given the results of DAMA [31] and Cogent [32] experiments. 

On the other edge of the mass range, the high mass regime, between 535 GeV and 50 TeV, 
is quite generic for scalar candidates of minimal dark matter models [25]. A detailed analysis of 
this mass range was given in [33]. They found that to satisfy the WMAP constraint, a small mass 
splittings between the components of the inert doublet is required, AM < 20 GeV. In this case, 
dark matter annihilates mainly into W + W~ , Z°Z° and hh. Previous analysis, including prospects 
for dark matter detection, were presented in Q3, 25,28 30] . A more recent study showed that the 
high mass regime of the inert doublet could be responsible for the DAMA signal through inelastic 
scattering [2"5] . 

The intermediate mass regime, M H o < Mw, is a particular feature of the inert doublet model 
which can not be reproduced in any other of the SU(2)l minimal dark matter scalar models. For 
these masses, dark matter annihilations are dominated by bb and by the three-body final state 
WW*, and the correct relic abundance can be obtained with or without coannihilations between 
H° and A (collider constraints on new charged particles prevent coannihilations with i? ± ). Given 
the richness and the simplicity of the physics involved, various studies have been devoted to the 
inert doublet model in this precise mass range. The first analysis of this parameter space were 
carried out in [T31[T4]. Further on, it was shown that the annihilation of inert Higgs particles 
could give rise to significant gamma ray lines [T5]. In [TlHT()ll2"8"ll3"0"] . the prospects for direct and 
the indirect detection of dark matter were considered. Constraints from LEPII data were derived 
in [T7], whereas prospects for detection at colliders were analyzed in [51H55] , Recently, it was 
demonstrated that dark matter annihilation into the three-body final state WW* plays a crucial 
role in this regime, modifying in a significant way the viable parameter space of the model and 
the prospects for its detection [20] . 

Notice, in particular, that none of these known viable regions includes the mass range Mw < 
M H o < 500 GeV. It is our goal to demonstrate that the lower part of that range, up to about 
160 GeV, is another viable region of the inert doublet model. Along this new viable region, 
cancellations between the different diagrams that contribute to dark matter annihilation into 
gauge bosons are required to satisfy the relic density constraint. 



4 Cancellations in the M H o > Mw regime 
4.1 Analytical Results 

Above the W threshold, M H o > Mw, dark matter annihilation into W + W~ becomes kinetimat- 
ically allowed, with the result that the total annihilation cross section tends to be rather large. If 
we estimate it from the pure gauge diagram -figured] (a)-, we get that 

(s -wavc) _ g 4 nrr 2 — (o , ^UMk ~ M w) \ n A , 

which, for M H a ~ 100 GeV, gives about 2 x 10 _24 cm 3 /s. And a similar number is found for 
the annihilation into Z°Z°. These annihilation cross sections are both much larger than those 
required to obtain the correct relic density (av ~3x 10 _26 cm 3 /s), so it seems difficult to satisfy 
the dark matter constraint in this mass range. In fact, previous analysis have claimed that it is 
not possible to do so. In [13], they stated that the inert Higgs could only provide a subdominant 
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(a) (b) (c) 

Figure 1: Feynman diagrams contributing to dark matter annihilation into W + W~ in the inert 
doublet model. The exchange diagram obtained from (c) is not shown. Analogous diagrams, with 
replaced by Z° and by A , contribute to the annihilation into Z°Z°. 



component of the dark matter, while in [14] they found the relic density to be too suppressed to 
be consistent with WMAP unless M H o > 600 GeV. 

Here, we show that those conclusions are incorrect. Taking advantage of the interference among 
the different diagrams that contribute to dark matter annihilation into gauge bosons, it is possible 
to satisfy the relic density constraint for M H o > My/- Thus, opening up a new viable region of 
the inert doublet model. 

Let us now explain how these interference effects come about. The diagrams that contribute 
to dark matter annihilation into W + W~ (or Z°Z°) are shown in figure [TJ Their amplitudes are 
given by 

iM p = i^ e *(pg) e *"(p 4 ) (15) 

= * s -m^L^ (p ^ m (16) 

iM t = ijA-^fps^M (17) 
iMu = irzW^CPs)*!^) (18) 

X 

where e*(p) is the polarization vector of the gauge boson with momentum p, pi and pi are 
the momenta of the H s, and s, t and u are the Mandclstam variables. In these equations, 
9v — g,g/ cos 9w and x = H^^A respectively for the annihilations into V = 14 7± , Z°. Typically, 
the first two contributions (the point-like interaction and the Higgs-mediated diagram) dominate 
the total annihilation cross section, so we will focus on those in our analytical discussion. Far 
from the Higgs resonance, their squared amplitude can be written as 




[( s -Mj) + 2\ L v^ 
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Since s ~ 4M H0 at low velocities, a cancellation between the two contributions occurs for 

A L w -2(M^ - (M h /2) 2 )/v 2 , (20) 

to which we refer as the cancellation condition. Hence, cancellations take place for Al > if 
M H o < Mh/2 and for Al < if M H o > Mh/2, For example, in the case of M H o = 120 GeV and 
M H o < M h /2, cancellations occur for Al ^ (0.04,0.5,1.2,2.0) when M h = (250,350,450,550) 
respectively, while for M H o > Mh/2, cancellations should occur for Al ~ [—0.35,-0.1] when 
Mh G [121,200]. Notice also, that the cancellation condition is the same for both final states, 
W+W~ and Z°Z°, indicating that both processes will be simultaneously suppressed. The idea, 
then, is that in regions of the parameter space where the cancellation condition is satisfied, the 
annihilation rate could be much smaller and, accordingly, a larger relic density will be found. 
Perhaps, large enough to be compatible with WMAP. 

The existence of cancellations in the annihilation cross section into gauge bosons had already 
been noticed in [17] and in [20] . But the possibility of using them to obtain viable models above 
Mw had not been considered before. 

Let us also emphasize that the t- and w-channels that we have neglected in the above analysis 
will certainly affect the annihilation cross section into W + W~ and Z°Z°. We will see that in 
order to comply with the relic density constraint it is necessary to suppress their contributions by 
considering large values of Ma and Mu± ■ 

Suppressing the annihilations into W + W~ and is not enough to ensure a small dark 

matter annihilation cross section. In addition, we must make sure that the other relevant anni- 
hilation final states, that is hh and tt, do not give a large contribution to it. If M H o > Mh, the 
total s-wave contribution to the HqHo — > hh annihilation process is given by 

r™0 _ A| / (Mj - 4Af4 _ 2M*v*\ L + 8M> 2 A £ ) 2 

aVhh -^M* hq V Mh ° Mh (M*-6M 2 M 2 HO +8M Ho y ■ (Zi) 

Using M H o > Mh and the Al of equation ([20]) . we obtain that cw^" wave ^ is typically of 0(lO~ 23 )cm 3 /s, 
giving rise to a very suppressed relic abundance. To be compatible with the observed dark mat- 
ter density, then, the mass of the dark matter particle must be smaller than the Higgs mass 
(Mho < Mh), so that the annihilation channel into hh is not kinematically allowed. This limita- 
tion of the parameter space is not that strong in the inert doublet model, for the Higgs mass can 
be quite large. A very stringent constraint comes instead from the annihilation into tt. The total 
s-wave contribution to the HqHq — > tt process is given by 

( s -wavc) _ N c Xjm 2 (M H o - mjf/ 2 

M ho tt (4M| - Ml) ' [ZZ > 

where m t is the top mass and N c = 3 is the color factor. Again we find that the Al satisfying the 
cancellation condition yields a very small relic abundance for M^o > m t . Therefore, viable models 
must feature M H a < m t . In |33j . for instance, it was assumed that M H o 3> mt, Mh, so the effect of 
cancellations was overlooked in their study. As a result, the minimal value of M H a consistent with 
WMAP was found to be 535 GeV. From our analysis of the mass range Mw < M H o < 500 GeV, 
we conclude that models satisfying the relic density constraint could be found only along the lower 
part of it, specifically for M H o < m t , Mh- 
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Figure 2: The dark matter annihilation cross section (at low velocity) into W + W~ and Z°Z° as 
a function of Al for two different values of Mgo: 95 GeV and 110 GeV. In this figure Mao = 
M H ± — 400 GeV and Mh = 130 GeV. Notice that cancellations take place simultaneously in both 
channels, W+W~ and Z a Z°. 



When evaluating numerically the annihilation cross sections and the relic density, we always 
take into account all the relevant processes, and all the diagrams that contribute to a given process. 
The analytical expressions derived in this section are not used for that purpose. We will verify, 
nevertheless, that the viable models always satisfy the cancellation condition we have obtained. 

These cancellation effects we have described are illustrated in figure [2j which shows the dark 
matter annihilation cross sections into W + W~ and Z°Z° as a function of \l- For this figure 
we set Mh — 130 GeV, M H ± = M^o — 400 GeV, and we consider two different values of M H o: 
95 GeV and 110 GeV. Equation (|20|) tell us that the cancellations should take place for for negative 
values of Al = —0.16 and —0.26 respectively which is in good agreement with the full numerical 
treatment of the annihilation cross section. We observe that the larger Mjjo the larger the value of 
| Ax, | required to obtain cancellations. Notice from the figure that cancellations may in fact reduce 
significantly the annihilation cross section and that this reduction takes place simultaneously in 
both channels, W+W~ and Z°Z Q . 

Since ilh 2 oc 1/crv, we expect these cancellations in the annihilation cross section to increase 
significantly the inert Higgs relic density, opening up the possibility of finding viable models for 
Mjfo > Mw- Next, we study, with several examples, the dependence of the relic density on 
the parameters of the models and show that there are indeed viable models in this region of the 
parameter space. 
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Figure 3: The relic density as a function of Xl for several values of M A a — M H ±: 200, 300, 400, 
500 GeV. In this figure M H o = 100 GeV and M h = 120 GeV. 



4.2 Examples 

Figure [3] shows the relic density as a function of Xl for M H a = 100 GeV, Mh = 120 GeV and 
different values of M^o and M H ± . For simplicity we consider models with Mao = M H ± and we 
let them vary between 200 GeV and 500 GeV. The horizontal band shows the region compatible 
with the observed dark matter density, 0.09 < Vlh 2 < 0.13. The effect of cancellations on the 
relic density is apparent in this figure. They increase tth 2 significantly as Xl approaches about 
—0.2, where the relic density reaches its maximum value. This value depends slightly on the value 
of M^o = M H ± but is in any case close to 0.1. Far from this maximum, say for Xl — —0.4 
or Xl = 0, the relic density is two orders of magnitude smaller. From this figure we also learn 
that cancellations are necessary but not sufficient to ensure a relic density compatible with the 
observations. In addition, the mass of the charged scalar (and of the CP-odd scalar) should be large 
enough to suppress the if ± (and ^4°) mediated contributions. If, for instance, M H ± < 200 GeV 
cancellations are still present but the relic density is always below the observed range. For the 
other values of M H ± we consider, 300, 400, 500 GeV, it is always possible to satisfy the dark matter 
constraint for a certain range of Xl- And the allowed values of Xl depend only slightly on the 
precise value of M H ± = M^o . This first example demonstrates that in the inert doublet model it 
is indeed possible to satisfy the relic density constraint for M H o > M\y- 

The dependence of the relic density on M H o and Mh is shown in figured] They both display 
the relic density as a function of Xl for M^o = M H ± = 400 GeV. In (a), Mho was set to 
100 GeV and four different values of M h (120, 130, 140, 150 GeV) were considered. As before, the 
horizontal band shows the region compatible with the observed dark matter density. The effect 
of cancellations are evident: they allow to satisfy the dark matter bound for Xl approximately 
between -0.15 (for M h = 150 GeV) and -0.22 (for M h = 120 GeV). As expected from equation 
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Figure 4: (a) The relic density as a function of A^ for different values of M%. For this figure 
we set M H o = 100 GcV and M A o = M H ± = 400 GeV. (b) The relic density as a function 
of Xl for different values of Mho. The other parameters were taken as Mh = 130 GeV and 
M A o = M H ± = 400 GeV. 



(|20l) , the value of Xl that gives the maximum relic density moves toward less negative values as Mh 
is increased. Because the Higgs-mediated contribution to dark matter annihilation is proportional 
to Xl, the Higgs-dependence of the relic density disappears for Xl = 0, as observed in the figure. 
At that point, the relic density is driven by the direct annihilation into gauge bosons (figure 
[2(a)) resulting in flh 2 ~ 10 -3 . We can conclude from this figure that for M H o = 100 GeV and 
M A o = M H ± = 400 GcV there is a significant range of Higgs masses for which the relic density 
can be consistent with the observations. 

In figure H (b) M h was set to 130 GeV and four different values of M H o (95, 100, 105, 110 GcV) 
were considered. We see that the position of the peak moves toward more negative values of 
Ai as M H o is increased, in agreement with the results found previously. To obtain the observed 
relic density, values of Xl approximately between —0.15 (for M H o = 95 GeV) and —0.27 (for 
M H o = 110 GeV) are required. These figures illustrate that the new viable region occupies a 
non-negligible volume in the parameter space of the inert doublet model, spanning a significant 
range in M H o, Xl, Mh, M A o and M H ±. 

So far we have illustrated the effect of cancellations on the relic density only for the light Higgs 
regime, Mh < 200 GeV. A remarkable feature of the inert doublet model is that it allows for a 
much heavier Higgs without contradicting electroweak precision data. As we will see, cancellations 
leading to viable models with Mho ^ Mw can also take place in the heavy Higgs regime of the 
model, that is for M h > 200 GeV. 

Figure [5] is analogous to figured] but for a heavy Higgs boson, Mh ~ 400 GeV. In both panels 
M H ± was set to 500 GeV and M A o to 450 GeV (they must now give a positive contribution to AT) . 
In (a) the relic density is shown as a function of Xl for M H o =120 GeV and four different values 
of Mh'- 360, 380, 400, 420 GeV. First of all, notice that in this case cancellations occur for positive 
rather than negative values of A^. As expected from the cancellation condition, the position of the 
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Figure 5: (a) The relic density as a function of Xl for several Higgs masses. In this figure we set 
M H o = 120 GeV, M A o = 450 GeV and M H ± = 500 GcV. (b) The relic density as a function 
of Xl for different values of Mjjo: 100, 110, 120, 130 GeV. The other parameters were taken as 
M h = 400 GeV, M A o = 450 GeV and M H ± = 500 GeV. 



peak moves toward larger values of Xl as Mh is increased. The values of Xl that give the correct 
relic density are between 0.5 (for M h = 360 GeV) and about 0.9 (for M h = 420 GeV). For Xl = 0, 
flh 2 becomes independent of Mh, with a value of about 6 x 10 -4 . In (b) the relic density is shown 
as a function of Xl for M h = 400 GeV and four different values of M H o: 100, 110, 120, 130 GeV. 
The viable region corresponds in this case to Xl between 0.7 (for M H o = 130 GeV) and 1.0 (for 
M H a = 100 GcV). We observe that the maximum value of fth 2 increases as M H o is decreased. 
In fact, if M H o =130 GeV the maximum relic density is as large as the WMAP lower limit. So, 
it is not possible to satisfy the relic density constraint (keeping the other parameters fixed) for 
M H o > 130 GeV. But if M H o = 100 GeV the relic density can be larger than the WMAP upper 
limit, reaching about 0.23. As a result, it is easier to satisfy the relic density constraint with a 
lighter inert Higgs. 

As we have seen in the previous figures, cancellations between the different diagrams that 
contribute to dark matter annihilation into gauge bosons indeed allow to satisfy the relic density 
constraint for inert Higgs masses above My/ (and below m t ). In other words, there certainly exist 
regions in the viable parameter space of the inert doublet model that had not been considered 
before. What we are lacking at this point is a more complete characterization of this new viable 
region. We would like to know, for instance, what is the largest dark matter mass (M H o) that is 
allowed, what are the typical values of Xl, Mh, M H ± and M A o, or what is the most promising 
strategy to detect the inert Higgs along this region. In the next sections we will address precisely 
these issues. 
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5 The new viable region 



So far, we have demonstrated the existence of a new viable region in the inert doublet model. 
This previously overlooked region features dark matter masses larger than My/ and is consistent 
with the relic density constraint thanks to cancellations between the Higgs mediated diagram and 
the direct annihilation diagrams into W + W~ and Z°Z°. In this section we perform a systematic 
analysis of this new viable region of the parameter space. To that end, we make use of Markov 
Chain Monte Carlo (MCMC) techniques to efficiently explore the full parameter space of the inert 
doublet model and find a large number of viable models in the region of interest to us. Then, we 
proceed to analyze such models in detail. 

5.1 The scan 

The inert doublet model contains 5 free parameters that can be taken to be Mjjo , M H ± , M40 , 
Al, and Mh- We allow these parameters to vary within the following ranges: 



80 GeV < Mjjo < 200 GeV (23) 

80 GeV < M H ± < 1000 GeV (24) 

80 GeV < M A o < 1000 GeV (25) 

-12 < \ L < 12 (26) 

114 GeV < M h < 600 GeV (27) 



MCMC tecniques allow us to scan freely and efficiently over these 5 parameters. Let us emphasize 
that these techniques have already been successfully used to scan and analyze the viable parameter 
space of other models of dark matter, see 38, 39] . Here we closely follow the procedure outlined 
in [38]. 

After scanning the above parameter space, we get a large sample of viable models, satisfying 
all the experimental and theoretical bounds and the relic density constraint. For our analysis 
we generated about 40000 such models within the new viable region. In the remaining part of 
the paper, we will work exclusively with this sample of models, analyzing its main features and 
studying its implications. 

5.2 Different projections of the viable region 

To facilitate the analysis and visualization of the 5-dimensional region where the viable models 
arc located, we will project it onto different two-dimensional planes. Figure El for instance, shows 
the new viable region of the inert doublet model in the plane (Mho To generate this kind 

of figures, we created a suitable grid in the two variables (M H o and Al in this case) and we 
place a point at a given value of them if there is at least one viable model in our sample with 
such coordinates. Sometimes we also differentiate (using symbols and colours) various ranges for 
other important parameters (the Higgs mass in figure [6|) that may provide additional information. 
Notice in figure [6] that the value of Myo over the new region extends from My/ up to about 
160 GeV. We argued in section 2] that M H o had to be smaller than m t — 171 .4 GeV but at that 
point we did not know how close to m t it could get. In Al, the range of possible values extends 
from —0.5 to about 3.0. A sharp contrast is observed between the models with Al < and those 
with Al > 0. In the former case (Al < 0), the maximum value of Mjjo is about 130 GeV and the 
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Figure 6: The new viable region projected on the plane (M#o,Al). The different symbols (and 
colors) distinguish among the possible ranges of Mh- The dashed line is the contour ^ = 0. 
Above that line /i| is negative whereas it is positive below it. 



Higgs boson is necessarily light, Mh < 200 GeV, as cancellations occur for M H o > Mh/2. In the 
latter case (Xl > 0), the maximum value of M H o can be larger, strongly depending on the range 
of Higgs masses. Thus, for 200 GeV < M h < 300 GeV the maximum M H o is only about 110 GeV 
whereas it can reach almost 160 GeV for the largest Higgs masses, Mh > 500 GeV. It is also 
apparent from the hgure, that for Al > the larger the Higgs mass the larger Xl must be. Given 
that cancellations are needed to obtain the correct relic density, it is not surprising that no viable 
models are found with Xl close to zero, as in that region the Higgs- mediated diagram becomes 
negligible. Finally, the dashed line shows the contour n\ = 0. From the figure we see that viable 
points with > feature A^ < and viceversa. 

The dependence of the new viable region on M H ± and M^o is illustrated in figure [7] Panel 
(a) shows the projection on the plane (Mho, M h ±) for different ranges of Mao. The trend is 
clear: the heavier the dark matter particle the larger M H ± must be. If M H a < 100 GeV, M H ± 
can take values within a wide range, but as M H o increases, higher and higher values of M H ± are 
required. If, for instance, M H o > 130 GeV, M H ± should be larger than about 400 GeV. And 
the viable models with the highest dark matter masses (M#o > 150 GeV) all feature a heavy 
(and a heavy A ), M H ± > 500 GeV. It is also observed in the figure that there is a correlation 
between the values of Mu ± and Mao . Large values of Mh± are associated with large values of 
Mj^o. This correlation is illustrated in more detail in panel (b) of the same figure, which shows 
the viable region in the plane (M^o,M H ±) for different ranges of the Higgs mass. In this plane, 
the viable models are concentrated along a band just above the M H ± = M^a line (the diagonal). 
In other words, most models feature M H ± > M^o but their mass splitting is never that large. 
Let us emphasize that this condition is not enforced by the dark matter constraint but rather by 
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Figure 7: (a) The new viable region projected on the plane (M H o ,M H ±). Different symbols (and 
colors) are used to distinguish the possible values of M^o. (b) The new viable region projected 
on the plane (Mao,M h ±). Three possible ranges of Mh are set apart with different symbols (and 
colors). 



electroweak precision data, see equations © and ([12"]) . We see then that both M^o and M H ± can 
be as large as 650 GeV and that the models with the highest dark matter masses (M H o > 140 GeV) 
require them to be larger than about 500 GeV. 

In figure [8] we show the new viable region in the plane (M^-o ,Mh). The Higgs resonance line, 
2M H o = Mh, is also shown for illustration. A priori we do not expect to find viable models close 
to that line, for resonant annihilations would drive the relic density to very small values. Points 
lying below the Higgs resonance line correspond to models with 2M^o > Mh and they all feature 
Xl < 0. From the figure we see that such points extend up to about 130 GeV in M H o and about 
160 GeV in Mh- They correspond, therefore, to what we call the light Higgs regime. Points above 
the Higgs resonance line (2M#o < Mh), in constrast, all feature a positive value of Xl- It can be 
observed that they include dark matter masses as high as 160 GeV and Higgs bosons heavier than 
200 GeV and extending up to 600 GeV. Hence, they lie within the heavy Higgs regime. Finally 
it is worth mentioning that no viable models are found for Higgs masses approximately between 
160 GeV and 200 GeV. 

We have in this way determined the position of the new viable region in the 5-dimensional 
space formed by M H o, -M40, M H ±, Xl and Mh- It is also useful to understand its position in the 
space of A parameters, which consists of A 3 , A4 and A5. The reason is that some constraints, such 
as the perturbativity and the vacuum stability bounds, are applied directly to these couplings, 
not to scalar masses or to A^. Since many models feature a rather heavy and ^ negative 
(equivalently Xl positive), we expect that a larger value of A3 will be required -see equation ([2]). 
That is what we observed in figure |9] It displays the new viable region in the plane (M#o ,A3) 
for different ranges of Mh- We see that A3 can indeed be large, with some models lying right at 
the perturbativity limit, A3 = 47r (dotted line). In the figure we observe that the models with 
the largest dark matter masses (M H o > 140 GeV), in particular, feature A 3 > 9. It must be 
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Figure 8: The new viable region projected on the plane (M H o,Mh). The line shows the position 
of the Higgs resonance region, 2M H a = M^. Different colors and symbols are used to differentiate 
models with Az, > from those with Al < 0. 



emphasized, however, that A3 is not that large over the entire new viable region. Models with 
A3 < 6, for instance, extend up to M H o = 120 GeV while those with A3 < 2 can reach up to 
105 GeV. Thus, if one wants to impose a more restrictive perturbativity bound on A3, this figure 
can be used to find out the resulting viable range of M#o . 

Figure [TU] shows the new viable regions in the planes (M h o,\a) and (Myo,As). A4 turns out 
to be negative in most models and reaches at most | A4 1 = 7, significantly below the perturbativity 
limit (not shown). The possible values of A4 move toward lower values as M H o is increased. 
Hence, for M H o = 100 GeV A 4 G (-1,-7) whereas for M H o = 130 GeV A 4 G (-3,-7). Notice 
also that A4 ~ — 5 for the models with the largest values of Mjjo. A5, which determines the mass 
splitting between H° and A , is by definition negative (M^o > Mho). We see, in panel (b), that 
its absolute value is never larger than 7 and only rarely exceeds 5, placing it comfortably within 
the perturbative regime. As with A4, the allowed values of A5 tend to move toward lower values 
(higher |As|) as Mho is increased. For the models with the highest dark matter masses, A5 is 
about —4.0. Consequently, a more restrictive pertubativity bound on A4 and A5, say A^s < 5, 
could be imposed without affecting the viable range of Mh° ■ 

We thus conclude this subsection regarding the projections of the viable region onto different 
two-dimensional planes. Before addressing the implications of the new viable region for the direct 
and the indirect detection of inert Higgs dark matter, we will discuss the issue of the fine-tuning 
required to obtain the correct relic density within this new viable region. Is it large, moderate or 
small? How does it compare against the fine-tuning found in other models of dark matter? 



16 



14 



Perturbativity limit - 




M o [GeV] 



Figure 9: The new viable region projected on the plane (M H o,X 3 ). The different symbols (and 
colors) distinguish among the possible ranges of Mh- The dotted line shows the perturbativity 
limit, A3 = 47r. 



5.3 Fine-tuning 

Given that in it cancellations are necessary to obtain the right relic density, one may think that 
this new viable region of the inert doublet model is disfavoured by fine-tuning arguments. Here, 
we quantify the fine-tuning of the viable models and show that that is not really the case. Many 
of them feature low to moderate fine-tuning and none of them has a very large one. 

The issue of the fine-tuning required to satisfy the dark matter constraint is not new. In [40] , 
where this fine-tuning was first quantified, it was proposed to use as a fine-tuning parameter the 
logarithmic sensitivity of the relic density to variations in the input parameters of the model. If 
di (i = 1, ...,n) are free parameters of a given dark matter model, the fine-tuning parameter (of 
the relic density) with respect to a,, a . is given by 

alogai 

And the total fine-tuning, total , is obtained by summing in quadrature the contributions of 

the different parameters of the model. Thus, in the inert doublet model we have that 

A n/l 2 itotal = ^A^ Vi2 Mh0 + A 2 nh2 MH± + A nh 2^ lA0 + A 2 nh2 Xl + A Qh 2 tMh . (29) 

Hence, the fine-tuning parameter is large if a small variation in the parameters of the models 
leads to a large modification of the relic density. If, for instance, Aq^ 2 , total < 10, a measurement 
of the parameters of the model at the 10% level will enable to compute f2/i 2 to within a factor 



17 



M H o [GeV] 



_J 

kid 



M H o [GeV] 



Figure 10: The new viable region projected on the planes {Mfjo^Xn) (left panel) and (M^o,As) 
(right panel). 



0(2). As a rule of thumb, one can say that the fine-tuning is small if tota i < 10, moderate 
for 10 < A n ^2 jtotal < 100 and large if Anha, total > 100. 

Figure [TT] shows the fine-tuning parameter as a function of M H o along the new viable region 
of the inert doublet model. Aq^2 total is typically smaller than 30, with plenty of models featuring 
A^ 2 .total < 10 over the entire range of M H o and no models with a fine-tuning parameter larger 
than 40. In particular, notice that there is no correlation between Anft^ total an d Mho. Low 
fine-tuning models, Ann 2 total S 10, exist for M^o just slightly above Mw as well as for Mh° > 
140 GeV. And the same holds for moderate fine-tuning models, A(ih 2 , total }t 20. 

These results can be compared with those obtained in previous works for one of the most studied 
models of dark matter, the Constrained Minimal Supersymmetric Standard Model (CMSSM or 
mSUGRA). It turns out that the viable regions of the CMSSM tend to be highly fine-tuned, with 
some regions, such as the focus-point and the funnel region, featuring a An/i= total greater than 
100 [40,31]. In contrast, in the new viable region of the inert doublet model we are studying, the 
fine-tuning parameter is never that high. 

6 Implications for detection 

In the previous section we described in detail the new viable region of the inert doublet model. 
Here, we examine its implications for the detectability of the inert doublet model via collider 
signals and the direct and the indirect detection of inert Higgs dark matter. 
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Figure 11: The fine-tuning parameter, tota i, as a function of M H o along the new viable region 
of the model. 



6.1 Collider signals 

The inert doublet model can be tested at colliders via the production and decay of inert particles 
or through the increase in the decay width of the Higgs boson [13]. Indeed, the decay processes 

h-> H H°,A A°,H + H- (30) 

may enhance the Standard Model Higgs decay width when Mh > 2Myo, with important impli- 
cations for Higgs searches at Tevatron and at the LHC. For the dark matter mass range we are 
studying, this increase is mainly due to the H°H° final state and it takes place only when Al > 
and Mh > 200 GeV, see figure [SI As a consequence, the Higgs decay width considered in the 
recent summary of Higgs searches at Tevatron below 200 GeV [42] is left untouched by the inert 
scalars for M#o > M\y- The LHC, on the other hand, is expected to perform direct measurements 
of the Higgs decay width with a precision of about 6% [43] for 300 GeV < M h < 700 GeV. 

Figure [12] shows the Higgs branching ratio into H°H° as a function of the Higgs mass for 
our sample of viable models. It lies approximately between 5% and 60%. In this figure, different 
ranges of M H o are distinguished. For models with M H o > 140 GeV, BR(/i — > H°H°) is observed 
to be about 10% or less. Branching ratios greater than 20% or so can also be found but only for 
models with the smallest dark matter masses, Mjjo < 100 GeV. In such case, the inert Higgs 
could be within the reach of LHC via measurements of the Higgs decay width. 

6.2 Indirect detection 

The indirect detection signatures of dark matter are determined, on the particle physics side, by 
the present annihilation rate and by the annihilation branching ratios. Figure [13] shows the present 



19 




M h [GeV] 

Figure 12: The decay branching ratio of the Higgs boson into the inert scalars as a function of 
Mh for different ranges of Mho. Only models with Ax > give a non-zero contribution to this 
branching. 



dark matter annihilation rate, cv, as a function of Mho along the new viable region (v is taken to 
be 10 -3 ). We see that most points are confined within a small interval between 0.5 x 10 _26 cm 3 /s 
and 4 x 10 _26 cm 3 /s. In other words, the present annihilation rate is always close to the typical 
value for a thermal relic, (cv) ~3x 10 _26 cm 3 /s. This outcome is a consequence of the fact that 
the annihilation of inert Higgs dark matter is dominated by its s-wave component, so crv remains 
unchanged between the early Universe and the present time. Besides, coannihilations, with A or 
H ± , play no important role in the region we are considering, as they would drive the relic density 
to even smaller values. 

In the new viable region, the inert Higgs annihilates mainly into bb, W + W~ and Z°Z° (the hh 
and tt final states are kinematically forbidden). Figure [T4l shows the annihilation branching ratio 
(BR) into these three final states as a function of M^o . We see, in the left panel, that the maximum 
value of BR(H°H° — > bb) decreases with M^o, going from about 0.85 for M#o ~ 80 GeV to about 
0.2 for M H a — 140 GeV. The minimum value of BR(H°H° -> bb), on the other hand, tends 
to remain constant at about 0.15. In particular, the bb final state is never negligible within the 
new viable region. The center panel shows the branching ratio into W + W~. As M H o crosses 
the Z° threshold, the annihilation into Z°Z° start playing a significant role and the branching 
into W + W~ is therefore reduced. It is also observed in figure [T4l that the range of variation of 
BR(H°H a -> W+W~) shrinks as M H o increases. For M H o > 120 GeV BR(H°H° -> W+W~) 
lies approximately between 0.3 and 0.5. Finally, the right panel shows BR(iJ°_ff° — > Z°Z°). It 
never exceeds 40%, reaching values between 0.25 and 0.35 for M H o > 120 GeV, and lower values 
for smaller M H o. Comparing the three panels we can deduce that bb, W + W~ and Z°Z° are 
indeed the dominant annihilation channels, as their combined branching is always close to 100%. 
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Figure 13: The present dark matter annihilation rate, <7v, as a function of Mho along the new 
viable region of the model. Notice that crv varies only over a small range, approximately between 
1 x 1CT 26 and 4 x lCT 26 cm 3 /s. 



The remaining annihilation channels, into light fermions, typically account for less than 10% of 
the annihilations. 

The present annihilation rate and the branching ratios that we have calculated determine the 
indirect detection signatures of inert Higgs dark matter. To illustrate the prospects for the indirect 
detection of dark matter along this new viable region, we compute the expected gamma ray flux 
from the annihilation of inert Higgs particles and compare it with present data from FERMI. 

The gamma ray flux from dark matter annihilations originating in a region AQ of the sky can 
be written as 

^=24^ rQ ^ JA ^' (31) 

where r Q = 8.5 kpc is the distance from the Sun to the Galactic Center, p Q = 0.39 GeV/cm 3 is 
the local density of dark matter [44] , dN/dE is the photon spectrum, and JAfi is given by 

J An = [ dfi(M) / - ( V . (32) 

J AQ Jlos r o V Po / 

For the galactic distribution of dark matter we consider a NFW profile [IS] , 

T ( V \ ^ 

PnfwM = Ps— H , (33) 

r V r s J 

where r s = 20 kpc and p s = 0.34GeV/cm 3 . For defmiteness we consider the observation region 
10° < |6| < 20°, for which there is available data from FERMI - see [IB] and the supplementary 
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Figure 14: The main annihilation branching ratios of inert Higgs dark matter as a function of 
M H o. The different panels show, from left to right, bb, W + W~ , and Z°Z°. 



online material. The photon spectrum from dark matter annihilations, dN/dE, depends on the 
annihilation branching ratios and has been obtained using micrOMEGAs. 

Figure [15] shows the differential gamma ray flux at 1 GeV (left panel) and at 10 GeV (right 
panel) for our sample of models. For E~ t — 10 GeV, it lies approximately between 5.0 x 10~ 6 
and 2.5 x 10~ 5 MeV/(cm 2 ssr), and at slightly smaller values (particularly at large masses) for 
1 GeV. Unfortunately, these fluxes are much smaller than those measured by FERMI [46], which 
are consistent with the predictions of cosmic ray models. It is thus not possible to obtain any 
constraints from this data. 

Let us stress that this somewhat disappointing result is consistent with previous analysis of 
FERMI data, see e.g. [47]. In them, it had been found that present data from FERMI can be used 
to rule out dark matter annihilation cross sections exceeding 10~ 23 or 10~ 24 cm 3 /s (depending on 
the dark matter mass and annihilation hnal state). Such cross sections are much larger than those 
expected in typical models of dark matter and, in particular, than those we have found along the 
new viable region of the inert doublet model. 

6.3 Direct detection 

Dark matter can also be detected via elastic scattering with terrestrial detectors, the so-called 
direct detection method. From a particle physics point of view, the quantity that determines 
the direct detection rate of a dark matter particle is the dark matter-nucleon scattering cross 
section. In the inert Higgs model, the H°N scattering process relevant for direct detection is 
Higgs- mediated, with a cross section, <th°n, given by 



where / ~ 0.3 is a form factor and m r is the reduced mass of the system. Thus, it is proportional 
to \\. Given the rather large values of Xl we found, see figure we foresee that cr H o N will be 




(34) 
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Figure 15: The expected gamma ray flux at 1 GeV (left) and 10 GeV (right) as a function of M H o 
for our set of viable models. The observation region is 10° < \b\ < 20° in both panels. 



significant over the entire new viable region. Figure \W\ shows the spin-independent dark matter- 
nuclcon cross section, ajjo N , as a function of M H a for our sample of models. <Jh°n li es m a 
narrow range between 3 x 10" 7 and 1CT 8 Pb. To ease the comparison with present experimental 
data, the bound from the CDMS experiment [37] is also shown (as a solid line) as well as the 
expected sensitivity of XenonlOO (dotted line) [48]. Points lying above the CDMS line are not 
consistent with present bounds whereas points above the XenonlOO line are within the reach of 
that experiment. As observed in the figure, a non- negligible fraction of the new viable region is 
ruled out by direct detection constraints. In particular, most models with Mh < 200 GeV (though 
not all of them) have a direct detection cross section exceeding the CDMS bound. No range of 
M H o is excluded, however. The good news is that all these viable models feature cross sections 
close to current bounds. They are all well within the sensitivity of XenonlOO. Clearly, direct 
detection is a promising way of probing this new viable region of the inert doublet model. 



7 Conclusions 

The inert doublet model is a simple and appealing extension of the Standard Model that can 
account for the dark matter. We have demonstrated, in this paper, the existence of a new, 
previously overlooked, viable region of this model featuring dark matter masses between Myy and 
about 160 GeV. In this region of the parameter space, the correct relic density can be obtained 
thanks to cancellations between the different diagrams that contribute to dark matter annihilation 
into gauge bosons. In the first part of the paper, we illustrated how these cancellations come about 
and analyzed the dependence of the annihilation cross section and the relic density on the different 
parameters of the model. Several examples of viable models were found both in the light Higgs 
regime, Mh < 200 GeV, and in the heavy Higgs regime, Mh > 200 GeV. In the second part of 
the paper, we made use of Markov Chain Monte Carlo techniques to obtain a large sample of 
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Figure 16: The spin-independent direct detection cross section as a function of Mfjo along the 
new viable region of the inert doublet model. Three possible ranges of Mh are distinguished by 
different symbols (and colors). The solid line shows the current bound from CDMS whereas the 
dotted line corresponds to the expected sensitivity of XenonlOO. 



models within the new viable region. These models were then analyzed in detail by projecting 
them onto several two-dimensional planes. From these projections we learned, for instance, that 
for Mh < 200 GeV the dark matter particle mass can be as large as 130 GeV whereas Xl is 
small and negative (0 > Xl > —0.5). Dark matter masses above 130 GeV (up to 160 GeV) can 
also be found but they feature heavy Higgs masses and positive values of Xl (3 > > 0). It 
was also noticed that the allowed values of M H ± and M^o increase with Mjjo along the new 
viable region, reaching Mao,M#± > 500 GeV for Mjjo > 140 GeV. In addition, we addressed 
the issue of the fine-tuning required to obtain the right relic density and showed that all these 
models feature a small or moderate value of the fine-tuning parameter. Finally, we considered 
the detection prospects of inert Higgs dark matter within the new viable region. For collider 
searches, the Higgs decay decay rate begins to differ from the Standard Model one for Higgs 
masses larger than 200 GeV. The extra contribution resulting from Higgs decays into inert scalars 
for Mh > 300 GeV might be observed directly at LHC. Regarding indirect detection, the present 
annihilation rate, the dominant branching fractions, and the expected gamma ray flux from dark 
matter annihilation were computed. But no bounds could be derived from current gamma ray 
data. The direct detection cross section was also calculated and found to be rather large. A 
fraction of the new viable parameter space, in fact, is already ruled out by present constraints 
from CDMS, and the remaining part will be entirely within the expected sensitivity of XenonlOO. 
Direct detection seems to be, therefore, the most promising way of probing this new viable region 
of the inert doublet model. 
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